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MULTI-SPEED SOLITARY WAVES 
OF NONLINEAR SCHRODINGER SYSTEMS: 
THEORETICAL AND NUMERICAL ANALYSIS 

FANNY DELEBECQUE, STEEAN LE COZ, AND RADA M. WEISHAUPL 


Abstract. We consider a system of coupled nonlinear Schrodinger equations 
in one space dimension. First, we prove the existence of multi-speed solitary 
waves, i.e solutions to the system with each component behaving at large times 
as a solitary wave. Then, we investigate numerically the interaction of two soli¬ 
tary waves supported each on one component. Among the possible outcomes, 
we find elastic and inelastic interactions, collision with mass extraction and 
reflexion. 
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1. Introduction 


We consider the following nonlinear Schrodinger system: 

J idtui + dxxUi + + P\u2\'^ui = 0 , 

1 idtU2 + dxxU2 + ^J.2\u2fu2 + /9|mi|^U2 = 0, 


(NLS) 


where for j = 1, 2 we have : R x R —>■ C, fij > 0, and /3 G R \ {0}. 

When /Ti = ^2 = /3, system (|NLSI) . also called Manakov system has been intro¬ 
duced by Manakov (see [21] for example) as an asymptotic model for the propaga¬ 
tion of electric fields in waveguides. In this particular case, it is to be noticed that 
the usual roles of x and t are inverted to study the evolution of the electrical held 
along the propagation axis. 

It has also been used later on to model the evolution of light in optical hber links. 
One of the main limiting effects of transmission in optical hber links is due to the 
polarization mode dispersion (PMD). It can be explained by the birefringence effect, 
i.e the fact that the electric held is a vector held and that the refraction index of 
the medium depends on the polarization state (see e.g mm- The evolution of two 
polarized modes of an electrical held in a birefringent optical hber link can indeed 
be modeled by (INLSI) in the case where = /i 2 and fj measures the strength 
of the cross phase modulation which depends of the hber (see ED)- Randomly 
varying birefringence is studied adding random coefficients in both nonlinearity 
and coupling terms of (|NLSIl (see for example [15] ) 

In higher dimensions, systems of nonlinear coupled Schrodinger equations ap¬ 
pears in various physical situations such as the modeling of the interaction of two 
Bose-Einstein condensates in different spin states. 

Systems of type (|NLSI) have also been studied from the mathematical point of 
view. When = ^2 = /?, in dimension 1, the system (|NLSp has the particularity to 
be completely integrable. Hence explicit calculations of solutions are possible and 
one can exhibit a variety of “truly” nonlinear solutions like solitons, multi-solitons 
or breathers (see e.g. the book m)- The integrability property is however not 
robust, and the slightest change in the parameters ^ 1 , /i 2 and j3 destroys it. Many 
works (see, among many others, [31 0 dH [29]) have been devoted to the study of 
the stationary version of (INLSp 

I dxx(t>i + -|-/3|02|^((>i = wi^i, 

1 dxxh + M2|</'2p(('2 + h = ^202, 


that one obtains when looking for standing waves solutions 


{ui,U2){t,x) = (e“^iVi(a;),e*‘^"V2(a:)). 


When standing waves exist, it is natural to study their stability and again many 
works have been devoted to this problem (see, again among many others, [THlISOllMl 
[28]). Existence and stability of standing waves are often proved using variational 
techniques. The analysis of 0 through variational techniques is very subtle and the 
introduction of new ideas is necessary to understand the full picture (see [H]). Note 
that (INLSp is Galilean-invariant. Hence a Galilean transform modihes a standing 
wave into a solitary wave traveling at some non-zero speed. 

Our goal in this paper is to provide a new point of view on the study of this 
system. We aim at understanding better the behavior in large times of solutions 
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starting at initial time as two scalar solitary waves carried by the two different com¬ 
ponents. We will use a mixture of theoretical and numerical tools, a combination 
seldom seen when dealing with this kind of problems. 

First, we propose to push further a study initiated in m on the multi-speed 
solitary waves of system (INLSI) and followed ud for a different nonlinearity in m- 
A multi-speed solitary wave is a solution of (INLSp which behaves at large time as 
two solitary waves. Here and as in we restrict ourselves to the case where the 
composing solitary waves are each carried on only one component of the system. 
In other words, taken independently, each component behaves as a scalar solitary 
wave at large time. Our first aim is to remove the high speed assumption under 
which the main result in m was proved. We therefore consider the system (INLSI) 
in dimension 1 and benefit from the fact that scalar solitary waves are in that case 
orbitally stable (see e.g [Tl]l. 

Our next aim is to investigate further the properties of multi-speed solitary 
waves solutions when they are crossing at positive time. Our theoretical result 
(Theorem[T|) indeed only guarantees existence of multi-speed solitary wave solutions 
to (INLSp when no interaction can occur at large time between the composing waves. 
But what happens when two solitary waves carried by different components collide? 
Due to the possible complexity of the phenomenon and the lack of appropriate 
theoretical tools to study it, we proceed the following numerical experiment. We 
take as initial data solitons on each components, both away from 0 but facing each 
other for the direction of propagation. Among the possible outcomes, we find elastic 
and inelastic interactions, interaction with mass exctraction, and reflexion. 


1.1. The theoretical result. Before stating our main theoretical result, let us 
give a few preliminaries. 

Let € H^{R) be the unique positive radial ground state solution to 

~ dxxQbj +ojQi^ — \Qu\'^Qu = 0 , Qui > 0 , Quj & ( 2 ) 

From simple calculations we note that the following scaling occurs: 

Quj{x) := yJTo Qi(v^a;). (3) 


For j = 1, 2, consider Wj > 0, 7 ^ £ M, Xj,Vj G M and define 




(4) 


The function Rj is a solitary wave solution to 

idtu + dxxU + = 0. (5) 

In this paper, we want to investigate the existence of solutions to (INLSI) where 
each component behaves like a solitary wave Rj solution to the scalar equation ([5|) . 
Our main theoretical result is the following. 


Theorem 1. Let > 0 and /3 £ R \ {0}. For j = 1,2, take Vj,Xj,jj £ R, 

ujj > 0 and consider the ground state profile Quij solution to ISD and the soliton Rj 
defined in ®. Then, there exist C > 0, Tq > 0 and {ui,U 2 ) solution to (INLSI) on 
the time interval [Tq, -I-oo) such that for all t £ [Toi +00), we have 




Wit)) 



^ Ce- 


H^xH^ 
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where a;» = 2;^ min{wi,a;2} and v^, = \vi — ^21- 

Remark 1. Compare to [m Theorem 1], the main differences are the following. 
Our result is valid for any speeds, whereas the one in m required a high speed 
assumption. We restrict ourselves to dimension 1 to have stable solitons (in m, 
any dimension was allowed). The overall proof strategy is similar, but in our case we 
need to perform several technical refinements which include in particular working 
with localized momenta and modulated waves. 

In addition, we are introducing the technical artefact consisting into introducing 
arbitrary constants in the definition (fTSl) of the global action. This is a new feature 
for this type of analysis, which is quite surprising as usually such a flexibility is not 
allowed by the algebra of the problem. 


The scheme of the proof is inspired by the one developed for the study of multi- 
solitons in scalar nonlinear Schrodinger equations in [nmniEsiES] ( see also [9] for 
a similar approach applied to Klein-Gordon equations). It consists in solving (INLSI) 
backward in time, taking as final data a couple of solitary waves (i?i (T”), i?2 (?"")), 
for an increasing sequence of times T" —+00. Thus we get a sequence (m",M 2) 
of solutions to (INLSI) on a time interval (—oo,T"'] such that (m"(T"'), U 2 (T”)) = 
(i?i(T"'),i?2(7’"))- We then have to prove the existence of a time Tq, independent 
of n such that for n large enough, is close to (i?i,i?2) on [To,T^]. The 

key tools at hand to prove Theorem [T] are 


• uniform in n estimates 


vte [To,r"], 


(uW) 




^ Ce- 




• a compactness argument that gives the existence of G iL^(]R) x 

iL^(R) such that converges strongly in iJ'*(R) (s € [0,1)) towards 


Remark 2. The method used to obtain Theorem[I]is a powerful tool to obtain sharp 
existence results for multi-solitons composed of ground states. Another approach 
relying on a fixed point argument has been developed for nonlinear Schrodinger 
equations in [nmg. This approach is very flexible and allows to prove existence 
of solutions more complicated than the ones in Theorem [T] like infinite trains of 
solitons or multi-kinks. The main drawback is that it always requires a large speed 
assumption. 


We also have tested numerically if the multi-speed solitary wave configuration 
was stable provided the starting waves are well-ordered and well-separated. In 
other words, we took the interaction to be small at the origin and the composing 
waves going away from each other. With such a well-prepared initial configuration, 
we remain close to a similar configuration in large time. This suggests that the 
multi-speed solitary waves are stable (no matter the coupling parameter). Note 
that this is expected due to the fact that each wave taken individually is stable. 
We have however no theoretical mean to verify this conjecture. Similar difficulties 
arise in the analysis of the stability for multi-solitons in nonlinear scalar Schrodinger 
equations (see e.g. [24]). 

1.2. The numerical experiments. We solve the system (INLS|) in one dimen¬ 
sion by adapting the time-splitting spectral method described in |7]. This method 
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is unconditionally stable, time reversible, of spectral-order accuracy in space and 
second-order accuracy in time, and it conserves the discrete total mass [5]. One 
can refer to [3] for other possible schemes and their properties. 

We will also compute the real valued ground state (minimizer of the energy 
on fixed mass constraints) of the system ([1]) using a normalized gradient flow 
approach. This will be used to make the comparison between the outcome of the 
interaction between two solitary waves and a solitary wave with profile 

As already mentioned, the experiment consists in taking as initial data the initial 
data of two solitary waves facing each other, each on one component. 

We considered four cases, the first one being the integrable case, where we expect 
the solitons after the interaction to move with the same velocity and amplitude. 
Apart when p,i = ^2 = /3, the system is not integrable, hence we do not expect 
pure elastic interaction between solitons. However, there are still regimes where 
we expect the outcome of interaction between solitons to be also a multi-speeds 
solitary wave, different from the input at two level: first, there are modihcations 
in the speeds and amplitudes of the composing solitons. Second, there is a loss of 
a bit of energy, mass and momentum into a small dispersive remainder. In certain 
cases, we have been able to identify the profiles of the outcome of the interactions 
as ground states of the stationary system O- 

The rest of this paper is organized as follows. In Section [21 we prove Theorem [1] 
assuming uniform estimates. In Section |3l we prove the uniform estimates. The 
numerical methods are described in Section IH and the numerical experiments are 
presented in Section |5l Appendix lAl contains the proof of a modulation result. 


2. Existence of multi-speed solitary waves 

This short section is devoted to the proof of Theorem [H assuming uniform esti¬ 
mates proved in the next section. 

In this section and in the next one, we assume that /ii, /i2 > 0 and j3 G \{0} are 
fixed constants and that we are given for j = 1,2 soliton parameters uij,Vj,Xj,'yj € 
K. Denote by Qui- and Rj the corresponding profile and soliton. 

We make the assumption that 


0 < III = —V2- (6) 

Since (INLSD is Galilean invariant, this assumption can be done without loss of 
generality. This will simplify calculations later on. 

Note that it follows from classical arguments (see [10]) that the Cauchy prob¬ 
lem for (|NLSI) is globally well-posed in the energy space x H^(R) and also 

in L^(K) X L^(]R). In particular, for any initial data (it®,!!®) G H^(R) x 
there exists a unique global solution (iti,it2) of (INLSI) in C(]R,7?^(]R) x ff^(R)) fl 
C^(R,H-^(R) X ff-^(R)). 

Let (T”) be an increasing sequence of times such that T" —>• -too as n —>• 
-too. Let (it",M2) be the sequence of solutions to (INLSI) defined by solving (INLSI) 
backward on (—oo,r"] with final data (m",M 2)(T”) = (i?i,i?2) (T”). The proof of 
Theorem [T] then relies on the following two ingredients. 

First, we have uniform estimates on the distance between the sequence (m",M 2) 
and the multi-speed solitary wave profile (i?i, R 2 ) ■ 
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Proposition 3 (Uniform estimates). There exist To > 0, no G N such that, for all 
n ^ no and for all t € [To,T"] we have 



(t) 



^ C , 


where, as in Theoreni[ll = 2 ^ niin{wi,^ 2 } and z;, = \vi — v^]- 


The proof of Proposition |3] is rather involved and we postpone it to Section [31 
The next ingredient is a compactness result on the initial data {uf,U 2 ^ (To). 
This result was already present in this form in m and we recall it without proof. 


Proposition 4 (Compactness). There exists G H^fR.) x such that, 

up to a subsequence, (u^jU^) (To) converges strongly towards (nJjU®) in H^fR) x 
i^®(R) for all s £ [0,1). 


With these two ingredients in hand, we can now conclude the proof of TheoremU 


Proof of TheoremlJl Let (ui, 1x2) be the solution on K of the Cauchy problem (INLSD 
with initial data (it^jU®) at t = To. By H^{R) x H^{R) boundedness and local 
well-posedness of the Cauchy problem in H^{R) x iJ®(]R) for all s € [0,1), we 
have weak convergence in H^{R) x H^{R) of {uf,U 2 ){t) towards {ui,U 2 ){t) for any 
t £ R. Combined with the uniform estimates of Proposition [31 this implies for all 
t £ [To, +00) that 




^ liminf 

mxm 




This concludes the proof of Theorem [TJ 


□ 


3. Uniform estimates 


This section is devoted to the proof of Proposition [3] In all this section, T” and 
(m", U 2 ) are given as in the beginning of Section [31 

3.1. The bootstrap argument. We first reduce the proof of Proposition [3] to the 
proof of the following bootstrap result. 


Proposition 5 (Bootstrap argument). There exist Tq > 0 and no G N such that 
for all n ^ Uq and for any tg £ [Tg, T"] the following property is satisfied. If for all 
t £ [tg, T”] we have 



it) 



(7) 


then for all t £ [tg, T”] we have 




( 8 ) 


The proof of Proposition |S] will occupy us for most of the rest of this section. 
We divided it into several steps. We first perform a geometrical decomposition of 
the sequence (uf,U 2 ) onto the manifold of multi-speed solitary waves in order to 
obtain orthogonality conditions. We then introduce an action-like functional, which 
turns out to be coercive due to our orthogonality conditions. This functional is not 
a conserved quantity, but since it is made with localized conservations laws it is 
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almost conserved. Using that property and a control on the geometrical modulation 
parameters, we are able to conclude the proof of Proposition O 

Before going on with the details of the proof of Proposition [Sj let us show how 
it implies Proposition [3] 


Proof of Proposition\^ Since we have (m”,U 2)(P") = (i?", i?2 )(’^") the final 
time T”, by continuity there exists a minimal time to such that for all t € [tojP"] 
we have 


We prove that to 
for all t G [to, T"] 


w- 


Ri 

R2 


(t) ^ (g) 

mxm 

= Tq by contradiction. Assume that to > Tq. By Proposition [Sj 
we have 



it) 



Therefore by continuity there exists too < to such that on [too,T”] estimate dH]) is 
satisfied. This however contradicts the minimality of to. Hence to = To and this 
concludes the proof. □ 


For the rest of Section [3l To > 0 and no G N will be large enough fixed numbers, 
and we assume the existence of to > To such that for all t G [to,T"’] the bootstrap 
assumption o is verified, i.e. we have 

1(3)«-(«:)«) 

Our final goal is now to prove that in fact dH) holds for all t G [to,T"]. 


^ e 


ffl X ffl 


( 10 ) 


3.2. Modulation. Let us start with a decomposition lemma for our sequence of 
approximated multi-speed solitary waves. 

Lemma 6 (Modulation). There exist C > 0 and functions 

ujj : [to,T"] — >■ (0,-|-oo), Xj : [to,T"] — >■ R, f/j : [to,T"] — >■ R, j = 1,2, 
such that if for j = 1,2 we denote by Rj the modulated wave 

R,{t,x) = - i,(t)), (11) 

then for all t G [to,T"] the functions defined by 



satisfy for j = 1, 2 and for all t G [to, T”] the orthogonality conditions 

(ej(t),.Rj(t))^ = (sjit),iRjit)'^_^ = (ej{t),d^Rj{t)^^ = 0. (12) 

Moreover, for all t G [to,T^], we have 


2 


E 


^|dtwj(t)|^ + \dtxj{t) 


da jit) + -j--ujj{t) 



( 13 ) 
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Remark 7. It is to be noticed that estimate (II3 clearly implies that, for Tq large 
enough, for j = 1,2, and for all t G [to, T"^], we have 

Xj{t) > > 2L and Wj(t) > 1152a;*. 

2v 2 

Moreover, a better estimate can be obtained for ojj and will be stated later on in 
Lemma [T^l 

This type of modulation result is classical in the literature dealing with solitary 
waves of nonlinear dispersive equations (see e.g. the fundamental paper of We¬ 
instein m for an early version or |24) for a recent approach). Its proof consists 
essentially in the application of the implicit function theorem combined with the use 
of the evolution equation to find equation (1131) for the evolution of the modulation 
parameters. We refer to the appendix for the details of the proof. 

3.3. Energy estimates and coercivity. In this subsection, we analyze the dif¬ 
ferent quantities that are conserved or almost-conserved in our coupled-vectorial 
problem. Remember that in the case of the scalar equation ([5]) the energy, mass 
and momentum, defined as follows, are conserved along the flow of ([5]): 

:= ^\\d^u\\l 2 - M{u) := ^||u||^2, P{u) := ]-Im f ud^udx. 

'J K 

The solution we chose of equation ([2]) is known to be the unique positive radial 
ground state of the action S := E(-, l)-|-a;M. Consequently, each soliton Rj defined 
by (HI) is a critical point of the scalar functional Sj defined by 

^3 ■= M + VjP. (14) 

Coercivity properties of linearizations of S),-like functionals are the key tool of the 

analysis of multi-solitons interaction (see for example [niEs]). _ 

In the vectorial case we are interested in here, the coupled system (INLSI) admits 
its own conservation laws. In particular, the mass of each component is preserved, as 
in the scalar case. However, the coupling does not preserve conservation of the scalar 
energy and momentum for each component and we only have conservation of the 
total energy (made of individual energies plus a coupling term) and total momentum 
(sum of the scalar momenta). More precisely, total energy, total momentum, and 
scalar masses of the whole system (defined as follows) are conserved quantities for 
the flow of system (INLSI) : 


£ (^2) E{ui, m) + E{u 2 , fj. 2 ) - ^ Jjui\'^\u 2 \^dx, 

V := P(ui) + P{U2), M, (^“ 1 ) := M(u,), j = 1, 2. 

In order to use the conservation of the momentum as in the scalar case, we here 
need to localize the momentum of each soliton, as was done in nans] for the scalar 
mass and momentum. Note that this was not needed for the analysis in m- Let 
us define the cut-off functions 

XLix) = x(J^), xi = l-XL, 







MULTI-SPEED SOLITARY WAVES OF NONLINEAR SCHRODINGER SYSTEMS 


9 


where L > 0 is arbitrary but fixed and x is a function such that 

0 ^ X ^ 1 on xi^) = 0 for a; < —1, xi^) = 1 for a; > 1, x^ ^ 0 on R, 
and satisfies, for some positive constant C and for all a; G R the estimates 
{x'{x)r ^ Cx{x), {x''{x)f ^ Cx'ix). 

Localized momenta are defined by: 

^loc (^,1 = ^ / i'^Aui + U2d^U2) Xl^x, j = 1, 2. 

\^ 2 / Z 

Remark that V = Note that since we are assuming (O the momenta 

above defined are localized around each composing solitary wave of the profile. The 
advantage of having made assumption ([5]) is that the cut-off does not depend on 
time. This will simplify our next calculations. 

In the sequel, we are interested in the following global action: 



= 8 



OJjit) 



+ Y. 

i=1.2 



+ Ci{l3, vi)M2 



+ C'2(/3, V2)Mi 



(15) 


where Cj{(3, Vj) are positive constants depending only on (3 and Vj and whose exact 
values will be decided later on. Note that the action implicitly depends on t via 
Cbj. It is to be noted that, in this work, we have the freedom to add these two 
coupled-mass terms that do not appear in the usual definition of iS-like functionals. 
This is a key point in our analysis. 

Let us now state in the following lemma several estimates related to the local¬ 
ization of Ri and R 2 and which will be of great use in the sequel. 


Lemma 8. For j = 1,2, if Tq is large enough, then for all t G [to,T"] and for all 
a; G R we have 

(\Rj{t,x)\ + \d,,Rj{t,x)\^ X^l\x) ^ C{1 + (16) 

n {\Rkit.x)\ + \dMt,x)\) < C(1 + |ui| + |u2|)2e-3V^’'*‘e-^l"l. (17) 

fc=l,2 

Lemma [5] follows from the support properties of the cut-off function and the 
exponential localization of the solitons profiles. Indeed, recall that in fact, the 
profile Q = Qi is explicitly known 

Q{x) = 2 sech(x), 

and it follows that Q and its derivatives are exponentially decaying, i.e. for any 
r] < I we have 

|Qxx| + |Q.| + |gKC'„e-’'l"l. 


Proof of Lemma O We prove only (I16|) , the proof of (I17|) following from similar 
(simpler) arguments. For simplicity in notation, assume j = 1, the case j = 2 being 
perfectly symmetric. Due to the exponential decay of the soliton profiles, we have 
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The cut-off function is supported on (—oo,i], and since for Tq large enough 
xi > 2L, for X S (—oo,L] we have (we recall here Remark[7]) 

\x-xi\^\x\, and |a; - ii| ^ 

As a consequence, we have 

In addition, as noticed in Remark [3 we have 

xi ^ — -=t^ and Cji ^ llb2u)^^ 

2v 2 

which implies 

|i?i(t,x)|xi(x) < C'e-3V^"**e-^l"l. 

The derivative dxRi is treated in the same way, with the only difference that |i;i| 
now appears in the estimate, due to the term in the definition of i?i. 

This finishes the proof. □ 

Lemma 9 (Expansion of the global action S). For all t € [to,T"] we have 


S 


Ui 

U2 

n 


= s 


Ri 

R 2 

Rfri 


U 


0(g-3yU7«.t)^ 


(18) 


coupled 


^2 


with 

Rfree 




1 = 1,2 


- I e'^Rjdx -I- Im I SjdxejXLdx 


1 


and 

"^coupled 


= C'i(/3,i;i)||e2|li2 -f C2(/3,U2)||£i||i2 


dx 


\e,f\R2f + \e2f\Ri\^) 

+-viRn [ £ 2 dxe 2 x\dx +-V 2 lm f £idx£ix\dx. 
^ IR ^ IR 


Proof of Lemmal^ First note that, according to the definition (fT^ of S, 


s {Z) = , 1 : 4(“>) 

j —1)2 


Kc 


- P{Ui) 


F C'i(/3,ui)M(ui) -I- C2{I3,V2)M{u2) - 


/3 


\ui\^\u2\^dx, (19) 


where Sj denote the same functional as Sj (see the definition (fTHl 'l with ujj instead 
of ujj. For j = 1,2, let us expand Uj{t) = Rj{t) + £j{t) in the expression of Sj. A 
simple computation leads to: 
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Now, as Rj is a critical point of the scalar functional Sj, we have 


S'{Rj) = 0, 


and thus 

where 


SjiRj+Sj) = S,iR,) + {S''iR,)e„e,)+0{\\e,r^,), 


( 20 ) 


1 


S”{Rj)ej,ej'^ = ^Wd^EjWl^ + ^ We^Wl^ + ^v^Xm / ejd^Sjdx 

dj. 


——TZe / e^R'^dx — / |£,p|-R, pdx. (21) 

2 Jr Jr 

Let us now develop the remaining terms in (1191) . As far as the momentum part is 
concerned, we write the expansion for j = I for simplicity: 


Ko 


Rl + £i 

R 2 + £2 


— p(i?i + £1) 


=-Xm / RidxRixlJx + Rn / R 2 dxR 2 x\dx 

Jr Jr 

- Imy ^.Ri5x£i + £i9x-Ri) x\dx+Xm J (^R 2 dxe 2 + e 2 dxR 2 ^ x\dx 

-Xm / £idxeix\dx +Xm / £ 2 dxe 2 x\dx. (22) 
Concerning the P coupling part in (fTg|l . we get: 

f I Rl H“ 1 ^ I R 2 “t“ ^2 P dx 
Jr 

= [ \Ri\^\R 2 \^dx^ 2 T^ [ (\Ri\^'^e 2 ^ \R 2 \^R~iei) dx 

+ J l-^2| H“ 1^2! |-^i| 4 dZe[EiRi)'Re{ 62 R 2 )^ dx 

+ 2 / (\£2\^R£{SiRi) \£i\‘^R£(^S 2R2)\ dx[ l^i p |£:2 (23) 

Jr^ ^ Jr 

Finally, the extra-masses terms in (1191) expand into 

M{R,+e,) = ^\\R^+e,\\l, = M(i?,)+£,)^ + M(£,) = PI{R,)+M{ej), (24) 

where we have used the orthogonality conditions (TT^ to obtain the last equality. 
In (I22I1 and (1231) , all terms containing a product of solitons or cut-off functions with 
different indices/exponents are of order 0(e-3v£^'«*4) by Lemma [H The terms 
containing a degree 3 or higher term in (£i,£2) are also of order 0(e-3V^''**) by 
the bootstrap assumption (ITOl) . Therefore, gathering (IT^ - (l20)) - (l^ -(l 
together gives: 


= 5 


Rl £1 
yi?2 + £2 
This concludes the proof. 


Rl 

R 2 


Rh, 


n 


coupled 


Ci(e 


— 3y/oJ7v*t\ 


□ 
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Lemma 10 (Coercivity of %). There exists A > 0 such that, for any to G [To,T"] 
and for all t G [fo, T"]; 

^2A(||ei||^, + ||e2||?,i), (25) 

'^coupled ^ A ^||ei||i2 + ||e2|lL2^ — A (^|| ||^2 + ||9a:£2|lL2^ , (26) 

and thus: 

^ (s) ^ ^ 



Proof. The proof of (1^51) is classical (see for example [HI US]) and we omit it. It 
remains to prove ( 1 ^ 51 ) . It is readily seen that for Cj{(3,Vj) large enough, we have 

C'l(/3,Ui)||£2||i2 +C2(/3,U2)||ei|li2 “ f (l^ll^l^ 2 |^ + |£ 2 H^i|^) dx 

^ 2 '*^1)11^211^2 + C'2(^,i;2)||ei|lL2) . 

As far as the momentum parts of Hcoupied are concerned, for example for j = I: 

^Tm / SidxSix^dx / \eidxeixWdx 

^ Jr ^ Jr 

^-yl|£i|l2l|5.^ill2 

^-^l|£i|li=-A||a.£ifL2. 

For C'(/3, Vj), j = 1,2 large enough, i.e such that 


we thus have: 

Hcoupied ^ A ^||ei||^2 + ||e2|li2^ — A ^||93;ei||^2 + ||(?xe2|li2^ • 
Combining estimates (1^ and (05)) on Hhee and Hcoupied finally gives (071) . □ 


3.4. Almost-conservation of the localized momentum. In this section, we 
investigate the conservation of the localized momentum and, inspired by [23], we 
state the following lemma: 


Lemma 11. There exists C > 0 (independent of L) such that if L and Tq are large 
enough then for all t G [to,T'^] and for j = 1, 2 we have 


^l(:m(o-^l ::)(T") 


L 


(28) 


Proof of Lemma [171 Let us prove the lemma in the case j = 1, the case j = 2 being 
perfectly similar. We first compute the time derivative of 

^7’ioc(^) = J dtUida:Ufx dx + ^Im J uidr,^tufx dx. (29) 
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Let US call term A and B respectively the first and the second term of the right 
hand side in equality 


A = ^Im J dtuAuix (0 dx, B = J uAdtUix dx. 

Using the fact that u{t) = ( ) is a solution to (INLSI) . we readily get (formally, 

\U2[t)J 

but this will be justified when integrating): 

dtUid^lT^ = idxxUidxUi + ifj,i\ui\^uidxui + il3\u2\'^uidxW (30) 

dxdtW = -idxxxW - - iPdx (\u 2 fwj ■ (31) 

About term A: Equation (1301) provides us with the following decomposition of term 
v4: 

A=l-Tie I dxxUidxUixiy) dx+^fiiBe [\uif uidxUiX (y) dx 

Z \L/ Z Jg \L/ 

+ ]-PTle [\u 2 fu 1 dxU 1 X (y) dx. 
Integrating by part in each term of A finally gives: 

= ■31/'®-“''“’'' (i) (D 


+ 


j\u2\^uidxuix(^^dx. (32) 


About term B: Equation (1311) provides us with the following decomposition of term 
B: 


B = - 


^Tie J uidxxxUiX dx-^^^iiTle. J uidoo ui'j x dx 

- 1 / 31 ^ J Uidx (^|u2|^Mr) X (0 dx. 

As for the A term, integrating by parts (several times if necessary) in each term 
finally leads to: 

® (l) ‘*" + 315(z) 

+ ^'“/j“‘U'(f)-i* 

+ 1/3"^ J \u 2 fu 1 dxU 1 X (0 dx + ^/3'^ J I“2 |^|mi|^x' (^) dx. (33) 

Combining (15^ and (1551) . and using the fact that x' and x"' are supported on 
[—L,L] lead to the following estimate: 

{\dxUi\^ + \ui\^ + dx+C J \u 2 \^ (^\uif + XLi^)dx. 




€ 


C 


l-L 
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Note that for x G [—L,L], Tq large enough, and using (fT^ we have (see Lemma [8] 
for similar arguments) 

\d,,Rj{t,x)+Rj{t,x)\ < ^ Ce-5\/^l“d ^ 

Note that C here depends on vj. Expanding now Uj = Rj + Sj, in (IMl) . using the 
above estimate and Lemma [H we obtain 


_d. 


rL{t) 


€ 


c 


0 {e 




H^xH^ 


The result follows integrating in time between t and T" and using the bootstrap 
assumption m- □ 


3.5. Control of the modulation parameters. We now give an estimate of the 
variations of uii with respect to time. This estimate is better than (1131) given by 
the modulation lemma. 

Lemma 12 (Variations of ojj{t)). For j = 1,2 and for all t G [toj T’"]? we have 

Proof. This estimate is due to the choice of the modulation orthogonality condition 
Rj{t)^ = 0, and the conservation of the mass of each component: 

0=\\n,\\l-\\R,\\l = \\R,\\l-\\R,\\l + \^^^^^^^ 

= {Cjj - Uj)— WQujWI + C>{\u}j - Wjp) + ||ej||2- 

OCJ 

Since < Oj ^^is concludes the proof. □ 

3.6. Conclusion. With the elements of the previous subsections in hand, we can 
now conclude the proof of Proposition [5] 


End of the proof of Proposition [31 Recall that we have made the bootstrap assump¬ 
tion (nni) and that our goal is to prove that for all t G [to,P”] we have in fact the 
better estimate 


Ri 

R2 


mxH^ 


1 

^ 2® 




(35) 


Let us first expand 



\R2{t) 

(um\ 


mxm 

\R2{t)J 


mxm 


(Mt)\ 

(Rim 

\Rm) 

\R2{t)) 





4- C ^ \ujj - ujj\+ ^ 0{\ujj - Uj\^). 
H^xH^ i=l,2 1=1,2 


By Lemma 1121 the part involving \ujj — ujj\ is controlled by the e-part. Hence to 
finish the proof it is sufficient to control e. Using successively the coercivity of Ti, 
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given in (|27jl . the expansion (|T^ of the global action S, the conservation of energy 
and mass, and the almost conservation of localized momenta ()28L we obtain 



2 

mxm 


< cn 





Wt)) 



L 


Therefore, choosing L large enough we obtain the required estimate (1351) . This 
concludes the proof of Proposition [5] □ 


4. Numerical schemes 

We describe here the numerical methods that we will be using in the next section. 


4.1. The time-splitting spectral method. We start by the time-splitting spec¬ 
tral method that we use to solve (INLSI) numerically. The equations are solved on a 
bounded interval I = {—a, a). We use a uniform spatial grid with mesh size h > 0 
and grid points Xk = xo + kh, k = 0,..., K, where iP -I- 1 G N is the (odd) number 
of grid points. Then h = 2a/K. The time grid is given by tn = n G No, 

where r > 0 is the time step size and to the initial time. We set {uj)/. := Uj(tn, Xk), 
where j = 1, 2, fc = 0,..., iG, and n G No. We split the system (INLS|) into the two 
subsystems 


idtUj = - /3\u3-j\'^Uj, j = 1,2, (36) 

idtUj — j — 


considered on and subject to some initial data. These subsystems are 

solved as follows. 

Step 1: Computing the evolution of (1361) we observe that the quantities \uj\'^ 
remain unchanged. Therefore, we “freeze” these values at time and solve 
the resulting linear ODEs exactly in the interval [!„,!„ -|- t/2], giving at 
time tn + t/2-. 


{ui)l = exp (i^(Mil(Mi)fer + /3|(u2)fep)) (ui)^, 
and analogously for (u2)J. 

Step 2: We solve (l37l) for j = 1,2 in the interval -I- r], discretized in 

space by the Fourier spectral method and solved exactly in time: 


iuj)T = 


1 


K + l 


K/2 

^ exp{-iTi'l/) {uj)^exp{wm{xk - Xo)), j = 1,2, 

n.=-K/2 


where Vm = 27rm /{xk — xq) and 


K 

(%)m = - ®0)) , 

1=0 



K 

T' 


Step 3: We solve (IMl) on -|- T/2,tn+i] using the discretization of Step 1 
with {uj)l* instead of {uj)^ and obtain (uj)'/'^^. 
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4.2. The normalized gradient flow. We will also need to compute the ground 
state solution of the following elliptic system with fixed masses 


= 0 , 


(38) 


j-dxx4>i + - P(i ^ - 

\ -dxx4>2 + W2(('2 - - P(t>l4>2 = 0. 

To that purpose, we use the normalized gradient flow. The problem can also be 
viewed as a nonlinear eigenvalue problem with u}i,u }2 being the eigenvalues, which 
can be computed from the corresponding eigenfunctions (j = 1,2): 


wj = 


JR 




■ Pfyl) dx 


(f)^jdx 


We solve 


by normalized gradient flow with given (01,02), such that: 

r r 

I 2 Q / 2 2 

/ (ji^dx = Oi and / 4 > 2 dx = 03. 

Js. Js. 


The standard gradient flow with discrete normalization consists in introducing an 
imaginary time in the nonlinear Schrodinger equations, thus looking at the imag¬ 
inary time propagation it —> —it) and after every step project the solutions such 
that the T^-norms are equal to (01,02). In the authors present the normalized 
gradient flow, prove it is energy diminishing, and propose numerical methods to dis¬ 
cretize it. Hereafter we adapted the normalized gradient flow for the given system 
and discretized it by a semi-implicit Backward Euler finite differences scheme. For 
i'Pj)k ~ 4’jit^j^k) being the discrete solution, Xk = xq + k ■ h the grid points with 
k = 0,1,... ,K — 1,K and the time sequence 0 < ti < t 2 < ■ ■ ■ < tn < tn+i < ... 
with T = tn+i — tn we have the following discretization of the normalized gradient 
flow: 


Step 1: We first solve on [tn,tn+i\ with initial data (j)j{tn,Xk)- 


i^j)l - i^j)k 


T 

with j = {1, 2}. As solution we get 
Step 2: is then normalized to get finally 






For t —>■ -1-00 we obtain the ground state solution {(j)i,(j) 2 ) of 
equal to (ai,a|) and frequency parameters 


with L^-norms 


5. Numerical experiments 

Our ansatz for the initial data in the next experiments is the following. 

= j = 1,2 (39) 

with defined in (|31). Without loss of generality, we may assume that 0 < iii = 
—V 2 (Galilean invariance), hence the soliton on the first component will be traveling 
to the right and the soliton on the second component will be traveling to the left. 
We will also assume that Xi = X 2 = 0 (invariance by translation in time and space) 
and to guarantee that our solitons on the first and second components are positioned 
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on the left and on the right respectively and at a sufficiently large distance, we will 
chose the initial time to be to = —10. 


5.1. Purely elastic interaction. The integrable case has been studied in depth in 
the book [1]. In particular, in that case the system is completely integrable via the 
inverse scattering transform and explicit solutions may be exhibited. For example, 
there exists a solution U = {ui,U 2 ) of (INLS|) which has the following behavior. At 
time t —>■ —c», 

uiit^x) - (v^(x - vit )), (40) 

U2{t,x) ^ (v^(x - V2t )), 

whereas at time t —>■ +c» the two components have interacted and the outcome are 
solitons of same speed and frequency but with a shift in phase and translation: 

ui(t, a;) - Mt, x) := (x - v^t - n)), (41) 

U2(t, x) ^ Mt, x) := (^ (x - V 2 t - T 2 )), 


where the shift parameters are given by the formulas 


Xi 


ln(|xi|) _ 
n =- t 2 = 

V^I 

■Cl - U2 + + v^) 

V 1 -V 2 + - y/uJ^) ’ 


ln(IX2|) _ Xi a _ 

_Vl- V2+ i‘2{y/uJl + y/iJ2) 
Vl - V2 - - y/w^Y 


Note that since the speeds and amplitudes of the solitons are not modified the 
interaction is considered elastic. Note also that we have a pure two-speeds solitary 
wave at both ends of the time line, without any appearance of dispersion despite 
the interaction. This is a characteristic feature of completely integrable systems. 

The parameters are chosen as follows: /ri = ^^2 = P = 1; = 5, a;2 = 1, 

Vi = 1, V 2 = —1, and Xi = X 2 = 0, thus we start with two solitons located at 
±10, respectively, which move to each other and observe them until tfmai = 10. 
Furthermore, the integration domain is / = (—20, 20) with K = 1024 spatial grid 
points and r = 10“^. 




Figure 1. Purely elastic interaction: Position densities 
|ui(t,x)p and \u 2 {t,x)p as functions of space and time for fii = 
/i2 = 1 and /3 = 1 
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We plot the solution ( 111 , 1 x 2 ) to (INLSp and observe that the two solitons remain 
unchanged after the interaction, excepting a shift in phase and translation. The 
shift in translation can be observed in Figure [T] Our numerical experiments are in 
good line with the theoretically predicted behavior. We have compared pointwise 
the numerical solution |ixip obtained by taking (l40l) as initial data at t = —10 and 
the theoretical outcome |(()ip given by (HD) . The results are shown in Figure [5] 



Figure 2. a) Plots of |ixi(10,x)p line) and |(^i(x)p (’.’ line), 
b) Difference of |iti(10,a:)p — \(j)i{x)\'^ 

5.2. Symmetric collision. In this experiment the outcome are still multi-speed 
solitary waves. After the collision we observe a small part of the soliton moving 
along with the soliton of the other component. We see that at the final time 
the solution fits well the ground state solution obtained by minimizing the system 
with fixed masses. Here we choose /xi = /X 2 = 1 and /? = 3. We have solved 
the system (INLSj) in the interval (—200, 200) with periodic boundary conditions, 
initial time to = ~10 and initial condition (1391) and the parameters wi = W 2 = 1, 
vi = 2, V 2 = —2, and xi = X 2 = 0, thus we start with two solitons located at ±10, 
respectively, which move to each other and observe them till tfinal = 40. Moreover, 
we use K = 4096 in space and time step r = 10“^. 

Here the interaction parameter /3 is greater than in the first experiment and after 
the interaction at time zero we see in each component a small part going with the 
other component, respectively (Figure [3]). 

At the final time t = 40 we look at the left side part of each component u~ = 
Uj ■ X[-oo,o]i and compute the T^-norms 

= f |it“(40, a;)pdx. 

JR 

Note that the mass of each component Mj{t) = ^ J^\uj{t,x)\'^dx is conserved, 
Mj(0) = Mj{t) for j = 1,2 and for all t G [—10,40]. At t = 40 we obtain (oi = 
V3^,a2 = yOTM). 

Finally we compute the ground state solution {cj)i{x),<j) 2 {x)) of the elliptic sys¬ 
tem (1331) by the gradient flow with constraints (3.893,0.069) and compare it to 
(ix)"(40,cc),ix^(40,x)). In Figure|4|we first shift <j)i{x + Xi) (with xi corresponding 
to max jitip) and (j) 2 {x+X 2 ) (with X 2 corresponding to max |iX 2 p) and then compare 
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Figure 3. Symmetric collision: Position densities |ui(t, x)p 
and \u 2 {t,x)\'^ as functions of space and time for /ii = /12 = 1 and 
/3 = 3 




Figure 4. Symmetric collision: a) Plots of |Mj"(40,a:)p (’- 
’ line) and + Xi)\^ line) b) Plots of |u^(40,a;)p and 

1 02 (a: + 2 : 2 )P, where Xj is the x corresponding to the max|Mjp, 
with j = 1 , 2 . 


the position density pointwise to the position density of the final solution located 
on the left of the origin (u)"(40,x),it^(40,x)) and realize that they fit very well. 

In this case, the interaction result into a new repartition of the mass and of the 
energy so as to approximate a ground state profile. 

5.3. Dispersive inelastic interaction. In this experiment we observe a loss of 
energy, mass and momentum in a small dispersive part located at the interaction 
place and a small dispersive part moving to the boundaries. Here we have solved 
the system (INLSI) with /3 = —1 and = /r 2 = 1 in the interval (—500,500) with 
8192 grid points for the spatial discretization, with periodic boundary conditions, 
initial time to = ~10 and initial condition (1391) with the following parameters 
vi = —V 2 = 2.7, = W 2 = 1 , and xi = X 2 = 0. 

Figure [S] shows the evolution of the waves. In FigurelSlwe see the position density 
of the initial soliton Mi(0,a:) and the one of the final soliton at t = 90, with lost 
mass and energy. 
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Figure 5. Dispersive inelastic interaction: Plots of \ui{t, x)\'^ 
(left) and \u 2 {t,x)\^ (right) as functions of space and time 




X X 

Figure 6. Dispersive inelastic interaction: Plots of 

|ui(0 ,x)p (left) and 11/1(90,x)p (right) 

5.4. Reflexion. In this experiment we observe a change of the sign of the speeds, 
thus the solitons are reflected after interaction. Here we have the same parameter 
as in the case of a dispersive interaction, /3 = — 1 and = /i 2 = 1, but the initial 
velocity is smaller, precisely we have the initial condition (I39p with the following 
parameters vi = —V 2 = 0.5, tui = 0/2 = 1, and xi = X 2 = 0. Since the velocity is not 
so large, and we observe the soliton until tfinal = 10, we solve (|NLSI) on (—20, 20) 
with 1024 spatial grid points, periodic boundary conditions and time step equal to 

T = 10-3. 

In Figure [3 we see the reflexion of the two solitons after interaction, thus only 
the velocities changed sign. In Figured we compare point wise the solution at the 
final time to the initial solution, and observe that it remains unchanged. 

These were various examples of the possible outcomes of multi-speeds solitary 
waves interactions. The above described situations may not be the only possibles 
and we do not claim comprehensiveness. 
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Figure 7. Reflexion: Plots of \ui{t,x)\'^ (left) and |u2(t,x)p 
(right) as functions of space and time 


|«i( —10, x) 1^, |«i(10, X-(-xi) 1^ 



Figure 8. Reflexion: Plots of \ui{—lQ,x — xi)\^ (’—’ line) and 
|ui(10,a:)p (’.’ line), where xi is the translation such that we can 
compare the two solitons 


Appendix A. Proof of the modulation lemma 

Proof of Lemma The proof of the modulation lemma is inspired from |23j and 
relies on the implicit function theorem. 

Let us define F : x K.+ x M x K. —> by 

/ (u-i?(w,7,i/),i?(a;,7,2/))2 \ 

F(u,a;,7,y)= (u - i?(w, 7, y), zi?(w, 7, i/))2 
\{u - R{uj, 7, y), d^R{uj, 7, y))^J 

where, for the sake of simplicity, we removed the i indexes and R{uj,"f^y) denotes 
the soliton i?(aj,7,y) = — y). Note that, 

V(w, 7, y) e M+ X K X M, F(i?(w, 7, y), w, 7, y) = 0 . 

The key idea here is to apply the implicit function theorem to function F at 
point (i?(a;, 7, y), w, 7, y) for fixed {u!,j,y) € M+ x R x M. Simple computations 
due to properties of Quj, mainly the fact that Qu: is radial and even, lead to 
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VF(i?(a;, 7, y),uj, 7, y) diagonal with diagonal terms a, b, c given by: 
a = —Tie I duiQuix - y)Qui{x - y)dx < 0 

M Jm 

b = - {dyR{uj, 7, y),iR{uj, 7, y))^ = —Re [ Quj{x - yfdx < 0 

d JR 

c = —\\dkQuj\\l 2 < 0. 

Implicit function theorem finally gives the existence of parameters Qj, 7^, Xj as 
functions of time. To prove that these functions are actually of class a stan¬ 
dard regularization argument is needed. We refer to [35] for more details on that 
argument. 

Now, in order to be more readable, we prove estimate m for j = 1. In that 
purpose, let us write the equation of evolution satisfied by ei, namely: 


idtSi + L{ei, 62 ) +M{ei,e2) 

= -{Qi - Y - i{vi - dtXi)dxQcj^. (42) 


where 


L{ei, £ 2 ) — dxx^i + + /3|.R2p^ £1 + 2y,iRe{Riei)Ri -f 2j3Re{R2£2)Ri 


and 

+ ^2/7,i7^(.Ri£‘i) + £^1- 

2 

It is to be noticed that the modulation terms cbi — ^— i9t7i, i9tWi, and Vi — dfXi 
appear in the right hand side of this evolution equation. The main idea of the 
proof is to take in (H3]) the scalar products of both sides of the equation with 
respectively Ri,iRi and d^Ri- We then use the orthogonality conditions (IT3]) as 
well as properties on . In the left hand side, we transfer the derivatives acting on 
£1 on the other side of the scalar product thanks to the modulation conditions (1131) 
(for the time derivatives) and integrations by parts (for the space derivatives). We 
finally use the equation at hand on i?i: 

idtRi + dxxRi + = 

V 4 J ^ fii 

(43) 

For more simplicity, we only develop the computations for the scalar product 
with i?i, the other cases are obtained using the same arguments. 
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Taking the scalar product with Ri in equation (j42)l leads to: 

^idtEi + L(ei,e2) + M{ei^e2),R\^^ ^ ll^ill 

-dtCjiTZ^ [ 

Jr 

- {vi - dtS:i)ne [ (44) 

Jr 

Right hand side of (1441) ; First, using equation dni) leads to: 


Tie / ie''^^^'"^''^~^'^^'>d^QcjiRidx = ——Im / duiQQi{x)Qcjj^{x - xi)dx = 0 
Jr vMi Jr 


and 


Re [ ie’‘^^'"^''"~^'^^'‘dxQcjiRidx = —^hn [ dxQi:ii{x)Qa,^{x - xi)dx = 0. 

Jr v^i Jr 

Thus, (imi reduces to: 

+ Li(ei,£ 2 ) +-^(ei,£ 2 ),^ = — ^ ~ 19*71^ l|RiII ■ (45) 


Left hand side of (gSl).- 

First, deriving modulation condition ^£i,Ri^ 


= 0 with respect to time gives: 


(^idtei.Ri^^ = (^ei.idtRi^ 

Let us now develop ^L(£i,£2),Ri 

^L(£i, £ 2 ), i?i^ ^ = ^£ 1 , dxxRi + 3/ri|i?i pi?i + /3|i?2pRi^ ^ + 2/3 ^£ 2 , |i?i ^ ■ 

Finally, ^A/'(£i, £2), Ri^ read: 

^A/'(£i, £ 2 ), Ri^ ^ + /3|£2pei,Ri^^ + ^(Miki P + /3|£2p)Ri, ^ 1 ) ^ 

+ ^^2/ii7^(Ri£i) + 2/37?£(R2£2)^ £ 1 ,^ 1)2 
Equation (|^ thus leads to the left hand side term: 

^£i,i9tRi + dxxRi + fill Rif Ri^ ^ + /9 |R 2 pRi^^ + ‘2‘fii ^£ 1 , |RipRi^^ 

+ 2/3 ^£2, |Ri|^R2^^ + ^A/'(£i, £2), ^1)2" 
Now, using equation (H 5 )) satisfied by Ri gives: 

^£i,i9tRi + dxxRi + Mi|Ri|^Ri)^ 

= (wi - ^ - dai) (£ 1 , Ri)^ + dtCji (^£ 1 , 

+ (ui - dtxi) (^£ 1 , 
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Finally, modulation condition 




= 0 leads to 


4 V 

+ (vi - dtx,) fsi, 

V VMi /2 

= —2^1 ^Ei, \ Ri\^Ri^ ^£i, —2/3 ^£ 2 , |i?l R2^ ^A/'(£i, £ 2 ), ^ ■ 

|46) 

Thanks to Lemma [ 8 l it is readily seen that most terms in this equation are of order 
0 (||£|| 2 )- Finally, (1461) can be re-written in a simpler way: 

2 

-Rill^ -b ai(t)) (wi - Y - 9ai) + a2{t)dtu}i + a3{t){vi - dtXi) = bi{t) 

where, for all t G [to,?"”], \ai{t)\ + \a 2 {t)\ + |a 3 (t)| -b |&i(t)| ^ C'||£(t)|| 2 . 

With the same kind of arguments, taking in (|42|) scalar product with iRi and 
idxRi respectively, we get two other equations that can be re-written as a linear 

/wi(t) - ^ -daiit)'' 

system solved by the ’’modulation vector” Mod(<) = dtCji{t) 

\ vi-dtxi(t).. 

This linear system takes the following form: 

(F(t) -b ^(t))Mod(t) = B{t) 

where 

/ pill 0 0 

F(t) = 0 j^du,Qcj^{x)Qcj^{x - xi{t))dx 0 

V f Pill 0 -ll^.g^Ji^ 

and, with the help of Lemma [51 for all t G [to,T^], P(t)||, p(t)|| ^ C'||£|| 2 , and 
hence (IT51) . □ 
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